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A dynamical transition separating intermittent and continuous flow is observed in a sandpile model, 
with scaling functions relating the transport behaviors between both regimes. The width of the active 
zone diverges with system size in the avalanche regime but becomes very narrow for continuous 



flow. The change of the mean slope, Az, on increasing the driving rate, r, obeys Az 



It has 



nontrivial scaling behavior in the continuous flow phase with an exponent 6 
only in terms of exponents characterizing the avalanches 6 = (1 + z — D)/(3 - 
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given, paradoxically, 
-D). 



Granular systems can exhibit continuous flow or inter- 
mittent avalanches depending on the driving rate. Exper- 
iments on rice piles have demonstrated that when grains 
are slowly added, transport of grains through the pile 
takes place in terms of avalanches of all sizes [jl[; these 
rice piles exhibit self-organized criticality (SOC) [§). 
Additional experiments have established that the trans- 
port process is dispersive, in the sense that the transit 
times of grains through the pile are broadly distributed 
||. A simple one dimensional "Oslo" model was pro- 
posed to mimic these experiments, and numerical simu- 
lations showed that the model exhibits SOC with disper- 
sive transport ||. This model was subsequently shown 
to represent a large universality class of avalanche phe- 
nomena including interface depinning, a slip-stick model 
for earthquakes Q , and maybe other sandpile models || . 

Although it is known that SOC can be reached only if 
the driving rate is very low, very little is known about 
the transition out of SOC as the driving rate is increased. 
In particular, for a given system size, there will be some 
driving rate at which the motion in the system never 
stops and the avalanches become infinite, signaling a new 
type of behavior which may or may not be related to the 
avalanche regime. Previously, Tang and Bak |j) mea- 
sured the increase in average height in the BTW sandpile 
model H as the driving rate was increased. They found a 
power law behavior (at high rates) for the average height 
vs. driving rate, but did not study the system at low rates 
(see in addition Ref. 0). Here we show, using the Oslo 
model, that there is an abrupt change in transport be- 
haviors distinguishing two regimes, an avalanche regime 
and a continuous flow phase, with a different critical ex- 
ponent for each one. This change is associated with a 
pronounced contraction in the width of the active zone of 
transport. We utilize the active zone behavior to relate 
the scaling coefficients in the continuous flow phase to 
the exponents characterizing the avalanches in the SOC 
state. 

The Oslo model is defined as follows: In a one dimen- 
sional system of size L, an integer variable h(x) gives the 



height of the pile at position x, and z(x) = h(x) — h(x+l) 
is the local slope. Grains are dropped at x = 1 with the 
opposite boundary open, i.e., h(L + l) = 0. At each time 
step, all sites are tested for stability. Each unstable site x 
with z(x) > z c (x) topples in parallel. In a toppling event 
at site x, h(x) —* h(x) — 1 and h(x + 1) — » h(x + 1) + 1. 
The key ingredient making this model different from pre- 
vious critical slope models [H is that the critical slopes 
z c (x) are dynamical variables chosen randomly to be 1 or 
2 every time a site topples. This randomness describes 
in a simple way the changes in the local slopes observed 
in the rice pile experiments || . 

We drive the model by adding grains at x = 1 at a uni- 
form rate r = 1/ AT, where AT is the number of lattice 
updates separating grain additions In the limit of 
zero driving, corresponding to SOC, no grains are added 
to the pile until the avalanche resulting from adding a 
sand grain ends and the system reaches a stable state 
with z(x) < z c (x) for all x. Clearly, for a finite system 
size L, well defined avalanches composed of activity sepa- 
rating intervals where the system returns to a metastable 
state are observed when the driving rate is nonvanishing 
but small. To be precise, for finite r an avalanche is 
defined as an interval separating two metastable configu- 
rations, where no topplings occur. The transition to con- 
tinuous flow can be viewed as a depinning transition away 
from the set of metastable states, where the avalanches 
become infinite. In general, if the time between addi- 
tions AT is much larger than the mean avalanche life- 
time in the SOC limit, (t(r — » 0, L)) the system will 
display intermittent avalanches; on the other hand, if 
AT <C (t(r — > 0,L)) the flow will never stop. Therefore, 
the relevant parameter for the avalanche to continuous 
flow transition is the ratio between these two quantities, 



R oc 



AT 



rL x 



where (t(r -> 0, L)) ~ L x with x = 1 + z - D ~ 0.19 for 
the boundary driven model. The dynamical exponent 
z ~ 1.42 ± 0.03 determines the cutoff in the lifetimes of 
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avalanches t co ~ L z and the exponent D ~ 2.23 ± 0.03 
determines the cutoff in the total number of topplings of 
avalanches s rn ~ L D in the SOC limit Elnt 



10000 




FIG. 1. Divergence of the mean size of the external 
avalanches with the rate, scaled with system size using 
x = 0.20 and D = 2.24. The fit parameters in Eq. (|l|) are 
a = 1.0, R c = 1.1, e = 1.2, and b = 1.64. The inset displays 
the same data but versus Ra j^ , with a double logarithmic 
vertical axis. 



In order to characterize the transition from intermit- 
tent avalanches to continuous flow we look at the mean 
size of the external avalanches (s(r, L)) out , which are the 
avalanches that drop grains outside the pile. Thus in 
determining the average, we only count avalanches that 
lead to an outflow. In the SOC limit the average exter- 
nal avalanche size resulting from a single grain addition 
is (s(r — > 0, L)) ou t ~ L D . As we increase the driving 
rate we find an extra, extremely fast divergence for these 
large avalanches, 

( A 

(s(r, L)) out ~ exp - 

\{r c - r) c 

when r approaches a transition rate r c . A collapse of the 
data is possible using the rescaled variable rL x and the 
scaling of (s) ou t in the SOC limit, which gives, 

(s(r,L)} ottt ~L D f(rL x ) . 

The consistency of this scaling ansatz is checked in Fig. 
|], allowing an accurate determination of the scaling ex- 
ponents, x = 0.205 ± 0.015 and D = 2.24 ± 0.02. Also, 
we show a power-law fit of the form 



log- 



<r,L)) c 
L D 



1 - 



lie 



- b 



(1) 



which locates the transition point at 
L 7 



with an exponent e ~ 1.2 ± 0.2. Note that it is the 
logarithm of the avalanche size which is diverging as a 
power law at the transition point. 

The key geometrical feature distinguishing the 
avalanche regime and continuous flow is the active zone, 
whose width is denoted as Al(t-) and is computed as 
X L (r) = ({h - {h)) 2 ) 1 ' 2 - Here h refers to the height of 
the pile. In the SOC state the active zone diverges with 
system size according to a power law Xl(t) ~ L x giving 
a rough surface for the pile, with x = D — 2; see Ref. 
0. One can argue that for small rates and fixed L the 
scaling of Al remains unchanged. In fact, we observe 
that the scaling of A^ remains unchanged over the entire 
avalanche regime. In the case of continuous flow, how- 
ever, the surface is smooth and the active zone is narrow 
with a width that is independent of system size. Thus 
at high driving, A^(r) — » A(r). Based on this fact, we 
propose the following finite size scaling ansatz, 

X L (r) ~ L*f{rL*), 

where the scaling function is constant for the avalanche 
regime and a decreasing power law for the continuous- 
flow phase, 



/(JZ)« 



constant for R < R r 



for R > R c , 



The exponent x/ x is obtained by imposing the inde- 
pendence of the active zone width on system size for 
fixed r > r c (L). Figure supports all of these results, 
giving a good estimation for the roughness exponent, 
X = 0.24 ± 0.01 and also xl x — 1-2 ± 0.1, in concor- 
dance with the values of x and X- Notice then that the 
properties of the continuous-flow regime depend on the 
exponents characterizing the SOC limit, without appar- 
ently introducing any new critical coefficients. 

These results suggest that the existence of an active 
zone that increases with system size is a unique feature 
of avalanche dynamics. In addition, the constant value 
of the roughness exponent x m the avalanche regime ex- 
plains the previously obtained constant value for the ex- 
ponent D observed in Fig. and shows the validity of 
the scaling relation D = \ + 2 for the entire avalanche 
regime. 

We also measured the dispersive transport properties 
of this system, by recording the transit time of each 
grain, T(r, L), which corresponds to the time a grain 
needs to travel through the entire pile. The distri- 
bution of transit times, P(T,L,r), is measured to be a 
(decreasing) power law for long times, with an exponent 
a ~ 2.2 ±0.1 that does not change with the driving rate. 
Therefore, the scaling relation a = D found in Ref. |ll]] 
seems to hold also for finite rate. A finite size scaling of 
this distribution is obtained scaling both axis with L z . 



r c (L) 



with R c ~ 1.1 ± 0.1 
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rL 0.20 

FIG. 2. Data collapse analysis for the active zone width. 
The exponents x = 0.20 and \ = 0-24 are used, and a 
power-law decrease with exponent x/x — 1.2 ±0.1 is obtained. 
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FIG. 3. Data collapse analysis of the scaling of the mean 
transit time with the driving rate, using x — 0.20 and 
z = 1.42. Power-law fits give exponents 1.00 ± 0.05 and 
a' = 2.2 ±0.1. 

It is possible to relate the average value of the transit 
time, (T(r,L)) to the width of the active zone. Using 
conservation of matter and incompressibility of the flow 
we get 



r(T{r,L)) oc LA L (r) 



(2) 



At small driving rate the average transit time of a grain 
scales as (T) ~ L D ~ X /r, in agreement with Ref. B. Con- 



sistent with this limit we can consider a general crossover 
function for both regimes, 

{T(r,L))~L z f(rL*) , 

using z = D - 1 + x with f(R) oal/R for R < R c . 

As a consequence of the behavior of Ai(r) above r c (L), 
the mean transit time scales linearly with L when r 
is kept fixed. Considering this limit in the crossover 
function the scaling function must also be a decreasing 
power law in the continuous flow regime, i.e., f(R) oc 
1 /R a for R > R c , with an exponent 

, z-1 



This behavior is demonstrated in Fig. |], with a' = 2.2 ± 
0.1, in agreement with the previous scaling relation. 

Clearly, as the system is driven harder the mean slope 
tends to increase. A power law behavior Az ~ r 1 / 9 may 
be observed, where Az is the increment of the mean 
slope with respect to the SOC limit, i.e., Az(r, L) = 
z(r, L) — z(r — > 0, L). We find that the observed power 
law varies depending on which regime the system oper- 
ates in. Since in the SOC limit, avalanches happen in- 
stantaneously, they do not contribute to the mean slope 
and then z(r — » 0, L) is the mean slope corresponding 
to the angle of repose of the pile. Therefore, for very 
small rates the mean slope will be given by z(r — > 0, L) 
plus the contribution of the avalanches. As the vertical 
scale of the avalanches is set by ~ L x , the mean dis- 
charge in slope during an avalanche will be proportional 
to Xl/L ~ 1/L 1 ^ X . For each addition AT this contribu- 
tion will have to be taken into account for a typical time 
(t(r -> 0,L)> and then 



Az(r,L) 



L AT L 1 -*-* ' 

which means that 9 = 1. We propose a scaling ansatz 

A*(r,L)~-^/(rL») , 

with f(R) oc R for small R. At r c the profile of the pile 
overcomes the maximum of the active zone fluctuations 
of the stable pile and then avalanches never stop. 

As the average transit time is linear with L for fixed r 
above the transition point, the grains see the same profile, 
independent of L, no matter how large the system size. 
Since z(r — > 0, L) is also independent of L, for large L, 
then Az is L-independent as well. Introducing this in the 
scaling ansatz we get our final result 



with 



f(R) oc R 1/e for R > R c 



x 1 ± z - D 



1-X 3-D 
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These results are displayed in Fig. |j, which determines 
6- 1 = 3.7 ± 0.15, that is, 9 = 0.27 ± 0.01, in agreement 
again with the obtained scaling relation. The fact of 
having 9 ^ 1 in the continuous flow regime can be used 
to argue that the infinite avalanche is not a superposition 
of finite avalanches, in contrast with the other regime. 
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FIG. 4. Data collapse analysis for the increase in the slope 
of the pile versus the driving rate, with 1 — \ — 0.76. The 
straight lines have exponents = 1.01 and # _1 = 3.67. 

It is worth mentioning that we have found a similar 
change in the value of 9 in the BTW sandpile model 0] 
and in the Manna model ]T^ | between low and high driv- 
ing rates. However, this transition is difficult to observe 
since in these cases the value of 9 in the continuous flow 
region is close to 1. This reconciles the results of Tang 
and Bak j|], valid for high rates, and a claim by Grass- 
berger and Manna in favor of 9 = 1 fDjfl . The transition 
we are reporting could also account for the different val- 
ues of exponents reported in the literature of interface 
depinning At present we are trying to extend our 

arguments to other models and compare the results with 
other approaches 

In conclusion, we have analyzed the transition from in- 
termittent avalanches to continuous flow that appears in 
a sandpile model when the driving rate is finite. For a 
finite but arbitrarily large system, the self-organized crit- 
ical state extends to finite driving rates with no change 
in the scaling properties associated with the roughness 
of the pile. Above a critical driving rate, the flow be- 
comes continuous. We argue that it is possible to relate 
the scaling properties of both regimes to each other with 
scaling arguments that are supported by results from nu- 
merical simulations. Remarkably, the nonlinear scaling 
properties of the continuous flow regime, where the pile 
appears physically to be completely different (having a 



narrow active zone) depend only on critical exponents 
associated with avalanches in the self-organized critical 
state. Since the Oslo model is a good approximate de- 
scription for transport and fluctuations of the profile of 
a real rice pile, it would be very interesting to test our 
theoretical predictions with experiments. 
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